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Abstract 



C/3 

^~>' A comprehensive modeling and simulation study of the first-order 

isotropic/smectic-A transition is presented and applied to phase diagram 
computation and two-dimensional spherulite growth. An approach based 
on nonlinear optimization, that incorporates experimental data (from 
12CB, dodecyl-cyanobiphenyl) , is used to determine physically realis- 
^ ' tic model parameters. These parameters are then used in conjunction 

VO , with an optimized phase diagram computation method. Additionally, a 

Cn ■ time-dependent formulation is presented and applied to the study of two- 

dimensional smectic-A spherulite growth. These results show the growth 
kinematics and defect dynamics of nanoscale smectic-A spherulite growth 
t*^ ' in an isotropic phase with an initially radial layer configuration. 

o 

oo 

O ■ 1 Introduction 



Liquid crystallinity and other forms of self-organization are key phenomena both 
technologically and in Nature. Liquid crystalline order ranges from liquids that 
jrt ' show some degree of orientational order to those that, in addition, show various 

types of translational order. The myriad of types of material that exhibit this 
behavior range from traditional low molecular mass molecules, currently used 
in display technology, to biological membranes composed of phospholipids [1] . 
To date, the main focus of liquid crystal research has been on the simplest of 
the class of mesophases where some degree of orientational order is present: 
the nematics. Research on mesophases that also show a degree of translational 
order, including smectics and columnar liquid crystals, has been less abundant. 



Recognizing the increasing importance of these mesophases, particularly in bi- 
ological systems, the need for practical methods to access the time and length 
scales at which these phenomena occur becomes more important. 

Experimental work in this general field has made great progress in the basic 
understanding of translationally ordered liquid crystals [2-6] . Nonetheless, it is 
currently infeasible to access much of the dynamic phenomena of translational 
phase-ordering processes. Recent experimental work has begun to address these 
issues [4, 5] , but theoretical approaches are currently the only way to access the 
length scales (nanometers) and time scales (nanoseconds) at which liquid crystal 
dynamics occur. The use of high-order models in conjunction with advanced 
numerical simulation techniques has shown a great deal of promise for theoretical 
study [7-10]. Recent computational advances have allowed for the possibility of 
simulation in greater detail than ever before. 

Utilizing a high-order Landau-de Gennes type model of the first-order isotropic/smectic- 
A mesophase transition [7, 11], the objectives of this work are: 



• 



to present a comprehensive approach to modeling and simulation of the 
first-order isotropic/smectic-A transition. 



• to determine phenomenological model parameters through incorporation 
of experimental data (from 12CB, dodecyl-cyanobiphenyl) . 

• to efficiently compute a phase diagram predicted by the model and pa- 
rameter set. 

• to study the two-dimensional growth kinematics and defect dynamics of 
an initially radial smectic-A sphcrulitc in an isotropic matrix. 

This approach builds upon previous work [12], which incorporates experimental 
data into the phenomenological model and derives equations for phase diagram 
computation. A time-dependent formulation [9, 10] and the nano-scale growth 
of an initially radial spherulite are presented. This work is organized as follows: 
a brief background on relevant types of liquid crystalline order is given (Section 
12. ip . the model and simulation approach are explained f Sections 12. 3112. 6p . and 
simulation results are presented and discussed (Section [3]) . 

2 Background and theory 

2.1 Liquid crystalline order 

Liquid crystalline phases or mesophases are materials which exhibit partial 
orientational and/or translational order. They are composed of anisotropic 
molecules which can be disc-like (discotic) or rod-like (calamitic) in shape. 
Thermotropic liquid crystals are typically pure-component compounds that ex- 
hibit mesophase ordering most greatly in response to temperature changes. 
Lyotropic liquid crystals are mixtures of mesogens (molecules which exhibit 



some form of liquid crystallinity) , possibly with a solvent, that most greatly ex- 
hibit mesophase behavior in response to concentration changes. Effects of pres- 
sure and external fields also influence mesophase behavior. This work focuses 
the study of calamitic thermotropic liquid crystals which exhibit a first-order 
mesophase transition. 

An unordered liquid, where there is neither orientational nor translational or- 
der (apart from an average intermolecular separation distance) of the molecules, 
is referred to as isotropic. Liquid crystalline order involves partial orienta- 
tional order (nematics) and, additionally, partial translational order (smectics 
and columnar mesophases). The simplest of the smectics is the smectic-A 
mesophase, which exhibits one-dimensional translational order in the direction 
of the preferred molecular orientational axis. It can be thought of as layers of 
two-dimensional fluids stacked upon each other. Schematic representation of 
these different types of ordering are shown in Figure [TJ 
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Figure 1: schematics of a) the isotropic, b) nematic, and c) smectic-A 
mesophases. 



2.2 Order parameters and the phenomenological model 

Theoretical characterization of mesophase order is accomplished using order pa- 
rameters that adequately capture the physics involved. These order parameters 
typically have an amplitude and phase associated with them. In order to char- 
acterize the partial orientational order of the nematic phase, a second order 
symmetric traceless tensor can be used [11]: 

Q = S (nn - ^l\ + l-P {mm - II) (1) 

where n/m/1 are the eigenvectors of Q-tensor, which characterize the average 
molecular orientational axes, and S/P are scalars which represent the extent to 
which the molecules conform to the average orientational axes [13-15]. Uniaxial 



order is characterized by S and n, which correspond to the maximum eigen- 
value(and its corresponding eigenvector) of the Q-tensor, S = |/j.„. Biaxial 
order is characterized by P and m/l, which correspond to the lesser eigenvalues 
and eigenvectors, P = —^ {fJ'm — l^i)- 

The smectic-A mesophase has one-dimensional translational order in addi- 
tion to the orientational order found in nematics. Characterizing this mesophase 
can be accomplished through the use of primary (orientational) and secondary 
(translational) order parameters together [16]. This is accomplished using the 
tensor order parameter ([T]) and the complex order parameter [11]: 

^ = Ve*"^ (2) 

where (j) is the phase, tp is the scalar amplitude of the density modulation. The 
density wave vector, which describes the average orientation of the smectic- 
A density modulation, is defined as a = V(/)/|V0|. The smectic scalar order 
parameter ip characterizes the magnitude of the density modulation, and is used 
in a dimensionless form in this work. In the smectic-A mesophase the preferred 
orientation of the wave vector is parallel to the average molecular orientational 
axis, n. 

A Landau-de Gennes type model for the first order isotropic/smcctic-A phase 
transition is used that was initially presented by Mukherjee, Pleiner, and Brand 
[7, 11] and later extended by adding nematic elastic terms [8, 17, 18]: 

/-/o - \a{Q:Q)^\b{Q-Q):Q + -^c{Q:Qf + ]^am'' + -^p\^>t 
-^^^ iQQ) \eQ ■■ {V^) (V**) 
+h,{VQ':VQ) + hi |V*|' + ^62 jV^^l' (3) 



A = ao(T - Tni) 
a = ao{T - Tai) 

where / is the free energy density, /o is the free energy density of the isotropic 
phase, terms 1-5 are the bulk contributions to the free energy, terms 6-7 are 
couplings of nematic and smectic order; both the bulk order and coupling of the 
nematic director and smectic density-wave vector, respectively. Terms 8-10 are 
the nematic and smectic elastic contributions to the free energy, respectively. 
The order parameters are defined in IBH]), T is temperature, T^j/Tai are the 
hypothetical second order transition temperatures for isotropic/nematic and 
isotropic/smectic-A mesophase transitions (refer to [19] for more detail), and 
the remaining constants are phenonicnological parameters. Further explanation 
and justification for the use of this high-order model can be found in [10]. 



2.3 Homogeneous free energy and parameter determina- 
tion 

Following past work [12], in order to compute the phase diagram from the free 
energy equation ([3]) a homogeneous uniaxially-ordered volume assumption can 
be used, resulting in a simplified point-volume free energy density: 

(4) 
where q = ^ is the magnitude of the wave vector and do is the equilibrium 
layer spacing. Note that a different definition of the Q-tensor ([T]) is used in this 
work than in ref. [12]. 

A major complication of applying this model is the determination of a suit- 
able set of model parameters. In order to overcome this challenge, a nonlinear 
programming formulation is derived which then allows to the application of 
nonlinear solution methods. Utilizing experimental data [19-21] and minimiza- 
tion criteria of the homogeneous free energy Q [7,22], a nonlinear optimization 
problem is formulated: 



obj == min {qf - q^) 



dfJTi) 
dX, 



^ 
> 
= (5) 



where the objective function obj minimizes the change in the layer spacing 
between the unknown bulk transition value, qt, and the known value at some 
minimum valid temperature for the model, qi. Each constraint is evaluated at 
the corresponding temperatures T;,, the bulk transition temperature, and Tj, 
the minimum valid temperature for the model parameter set. 

2.4 Phase diagram determination 

Utilizing the homogeneous free energy ([4]) and model parameters determined 
from the nonlinear optimization solution in Section 12.31 a phase diagram for 
the system can be computed. Figure [H shows a schematic of the phase diagram 
of a first-order isotropic/smectic-A phase transition. Coexistence regions are 
present, due to the first-order nature of the transition, where both the isotropic 
and smectic-A phase are either stable or metastable. 

Minimization of the homogeneous free energy ^ is computationally chal- 
lenging due to the presence of three degrees of freedom. This can be simplified 
by parameterization of the free energy equation [7] using minimization invariants 
and the assumption that the smectic-A phase exists. The resulting free energy 




Tai Tb 

Temoerature 



Figure 2: Phase diagram schematic of the first-order isotropic/smcctic-A phase 
transition where Tsh corresponds to the upper stability hmit of a super-heated 
smectic domain, Th refers to the bulk transition temperature (the free energy of 
the isotropic and smectic- A phases are equal), and Tai refers to the lower sta- 
bility limit of a super-cooled isotropic domain, also referred to as the theoretical 
second-order transition temperature [19]; the smectic-A phase is stable and the 
isotropic phase is unstable in region A, the smectic-A phase is stable and the 
isotropic phase is meta-stable in region B, the smectic-A phase is metastable 
and the isotropic phase stable in region C, and the smectic-A phase unstable 
and isotropic phase stable in region D. 



equation, parameterized as a function of the nematic scalar order parameter 
(see (dl), is: 

fA = AS\ + BSI + CSI + DSa + E 

A ^ - (lee"* + 96b2Se^ + lUbjS^ - lUbjfSc) (l2966^/3) "^ 

B = - (-966ie3 - 288bib2Se + 96bbj(3) (12966^/3)"^ 

C = - {{216bl - U4:ab2) e^ + {2mblb2 - ^32abl) S - A32ablf3) {l296blP)~^ 

D ^ - {A32abib2 - 2166?) e (12966^/3) "^ 

£; = - (324a2^2 _ 3240,^2^2 + 816^) (12966^/3)"^ (6) 

where Sa is the nematic scalar order parameter with the assumption that the 
smectic-A phase is stable/meta-stable. The minima of ^ are easily found in 
that they are the roots of a polynomial, -^^ . Once the minima are determined, 
the validity of the assumption that the smectic-A phase exists (at the specified 
temperature) can be tested by i) verifying that the computed Sa and the corre- 
sponding values of ip/q are real, and ii) verification of the minimization criteria 
of the full homogeneous free energy ([4]). The equations to determine q and ip 
from Sa are derived from the application of the minimization criteria to the 
homogeneous free energy equation Q : 

y 2eSA — 36i 

q = 



362 

6/3 



2 _ 4(55*;^ + AeSAQ^ - 3629^ - 66ig^ ~ 6a 



2.5 Time-dependent formulation 

The Landau-Ginzburg time-dependent formulation [23] is used to capture the 
kinetics of the phase transition. Due to the higher order derivative term in the 
free energy functional, a higher order functional derivative must be used. Addi- 
tionally, in order to utilize standard numerical solution techniques, the complex 
order parameter ^ is separated into its real and imaginary contributions: 

"$ = A + Bi (8) 

The general form of the time-dependent formulation is as follows [23]: 



dQ 



dt 



717 I I -If I (9) 




F = JdV (10) 



'V 



where ^ir/^J■s is the rotational/smcctic viscosity, / is the heterogeneous free 
energy density ([3]), and V the volume. As previously mentioned, a higher order 



functional derivative must be used due to the second-derivative term in the free 
energy Q: 



SF df d df 



df 



de dx, \d-l§-^ 



d d 



dx, dXj \ djr^ 



(11) 



where 6 corresponds to the order parameter. 

The use of the full Q-tensor in this time-dependent model does not neglect 
biaxiality as was assumed in Sections 12.3112.41 Substituting ([5]), the free energy 
([3]), and high order functional derivative piD into the time-dependent formula- 
tion ([9]) yields the closed set of simulated equations: 



dQ 

dt 

dt 
dt 



a*Q -b*{Q- qf^ + c*{Q:Q)Q~6* {A^ + B^) Q - ^e* (VylV.4 + VBVBf^ 



- {a* A + /3* (^2 + ^2^ ^ _ j*^ (g : g)] + V 

- {a*B + /3* (A^B^) B - 5*B {Q : Q)] + V 



blVA - e*Q ■ VA - -b;V (V^A) 



blVB - e*Q ■ VB - -6* V (V^B) 



where the asterisk denotes an nondimensionalized value, the superscript ST 
denotes the symmetric/traceless portion of a tensor, and /i* is the ratio of the 
smectic and rotational viscosities. The nondimensionalized model parameters 
are as follows: 



bl^ 



ao 
bi 
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(13) 



where I is the simulation-specific imposed length scale. 



(12) 



S* ^ 



aoAT 



2.6 Simulation conditions and computational approach 

A square computational domain with imposed length scales of / = 9.75 x 10~^ fim 
(approximately 25 layers) and I = 1.95 x 10~^/^to (approximately 50 layers) were 
used in two separate simulations. Referring to Figure [3^, both symmetry and 
Neumann boundary conditions were used to simulate bulk conditions. Symme- 
try conditions for the Q-tensor ([l]) must take into account vector symmetry, 
which results in the following boundary conditions [24] : 



dQ 



dx 



XX p,_ ^^yy 



dy 



— U, L^xy — ^yx U 



(14) 



The initial condition for both simulations was a smectic-A spherulite in an 
initially radial layer configuration (see Figure [31d) . The radius of the spherulite 
was initially set to vq — 12.0nm, a value on the order of the layer spacing at 
330-ftr (Section mH) , do = 3.9nm [20]. The initial value used for S, ip, and the 
layer spacing correspond to the homogeneous values at T = 330K, determined 
from the computed phase diagram. The Heaviside step function was used to 
generate the initial spherulite. The constraint that the spherulite does not 
impinge on the domain boundaries was verified post-simulation. 

A commercial package, Comsol Multiphysics, was used to solve the time- 
dependent model P^ . Quadratic Lagrange basis functions were used for the 
Q-tensor variables and quartic Hermite basis functions used for the complex 
order parameter components. Standard numerical techniques were utilized to 
ensure convergence and stability of the solution. This platform does not support 
adaptive mesh refinement, thus a uniform mesh was used with a density of 
approximately 14.8 nodes/nm^ for the 25 layer simulation and 3.7 nodes/nrn? 
the 50 layer simulation. Previous simulations using this model and numerical 
method have shown good agreement with both past experimental and theoretical 
findings [9, 10]. Additionally, exhaustive past work using this numerical method 
and the Landau-de Gennes model for the first-order isotropic/nematic phase 
transition [25-28] has served to further validate this simulation approach. 

3 Results and discussion 

3.1 Phase diagram computation 

The algorithm for phase diagram computation presented in Section 12.41 was im- 
plemented using a high-level numerical programming language [29] . The result- 
ing phase diagram is shown in Figured! refer to the phase diagram schematic. 
Figured for a detailed explanation of the features. 

The experimentally observed bulk transition temperature and nematic scalar 
order parameter [19] are exactly matched using the parameter determination 
method presented in Section 12.31 Additionally, the smectic layer spacing on 
the order of 3.9nm [20] is also well reproduced below the lower stability limit 
of the isotropic phase, Tai- Above this temperature the model predicts a layer 
spacing trend that increases approaching the super- heated stability limit of the 
smectic-A phase. 

3.2 Two-dimensional spherulite growth 

Simulation results from the growth of an initially radial spherulite are shown 
in Figure [5] from the 50 layer simulation (see Section 12. 6p . Past work on the 
isotropic/nematic mesophase transition [14, 24-28, 30] and smectic-A filament 
growth [31] provide a great deal of insight into the two general growth processes 
observed: shape dynamics [15] and self-similar growth [32]. As will be shown, 
these simulation results definitely show that the growth of an initially radial 
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Figure 3: a) (left) a schematic of the simulation domain showing boundary 
conditions and spherulite placement b) (right) Schematic of the initially radial 
spherulite layer configuration used as an initial condition for simulation. The 
center of the radial spherulite is initially assumed isotropic. 
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Figure 4: The computed phase diagram for the 12CB-based parameters. Note 
that q, the wave vector magnitude (see (HJ), is scaled by 3 gxio-iim I^*^]- "^^^ 
experimental bulk transition temperature for 12CB is 331.35 and the value of 
the nematic scalar order parameter at that temperature approximately 0.61 [19]; 
The model parameters used in this computation and droplet growth simulations 
are as follows: Tni = 322.85^, Tai = 330.5i4:, ao = 2 x 10^;;;;^, b = 2.823 x 
10^;^, c = 1.972 X 10^;^, ao = 1.903 x 10^^^, p = 3.956 x 10^^, S = 
9.792 X 10^;^, e = 1.938 x IQ-^^pN, h = 1 x IQ-^^pN, 61 = 1 x IQ-^^pN, 
62 = 3.334 x lO~^'^Jni, and the ratio of the rotational and diffusional viscosities 
used was -^^ = 25. The rotational viscosity value um = 8.4 x 10"^-^^-^ was 
used post-simulation for estimation of the time scale [25] . 
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textured smectic-A spherulite in an isotropic phase follows a similar topological 
growth process as observed in its nematic counterpart [25-28]. 




Figure 5: Results from the 50 layer simulation (domain length scale, utilizing 
symmetry, is I = 3.9 x 10~^fim) of the growth of an initially radial textured 
spherulite where the surface corresponds to A = i?e(*); (i) t = 3A0^s, the 
self-selected formation of a +1 disclination occurs from the initially isotropic 
core (see Section [2^ (ii) t — 17.1/is, spherulite growth continues with the same 
core morphology (iii) t — 28.7^s, growth continues with the splitting of the 
+ 1 disclination into two -1-^ smectic disclinations and the creation of a single 
disoriented smectic layer in the core; layer configuration results from the 25 
layer simulation are indistinguishable from the 50 layer simulation, up to the 
maximum radius achieved with the 25 simulation, and are excluded for brevity. 



3.2.1 Shape dynamic growth 

The shape dynamic growth regime is dominated by transient texturing and 
interfacial forces. When the spherulite radius is on the order of the smectic 
coherence length: 



A 



bi 



(15) 



long-range energy effects (gradients in molecular and layer orientation) are dom- 
inant over short-range energy effects (gradients in the bulk order parameters S 
and tjj). As the spherulite radius increases, interfacial anchoring affects the 
core of the spherulite less and short-range energy gradients become dominant 
in this region. This results in increased texturing and defect dynamics [13, 14]. 
Referring to Figure [SI this process is observed for an initially radial smectic-A 
spherulite. 

The shape dynamic growth regime of this spherulite morphology has been 
found to have topologically similar dynamics compared to the growth of ini- 
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Figure 6: Enlarged results of the spherulite core from the 50 layer simulation 
(figure length scale is / = idnm) corresponding to the full domain in Figure [5] 
where the surface corresponds to A = Re{'^); (i) t = 3.40/xs, the circular layer 
configuration in the vicinity of the +1 smectic disclination (ii) t — 17.1/xs, the 
splitting of the +1 disclination with the formation of a new layer (iii) t — 28.7/is, 
the pair of +^ disclinations are fully formed which results in a single disoriented 
layer undergoing expansion with a wave vector orthogonal to the (vertical) axis 
formed by the two disclinations. 



tially radial nematic spherulites [25-28]. In agreement similar nematic growth 
simulations [25-28], the initially imposed isotropic core is observed to undergo 
a self-selected transition to a +1 disclination, shown in Figure [7^. This ini- 
tial morphology minimizes elastic long-range energy through the formation of 
a high-energy defect. As the spherulite radius exceeds the smectic coherence 
length p5|) . long-range elastic gradients become energetically favorable com- 
pared to gradients in short-range order. Subsequently, the -1-1 disclination splits 
into a pair of -|-^ disclinations [14,33]. This is a topologically equivalent texture 
that minimizes short-range energy, where the energy of a defect is proportional 
to its the square of its strength s^ [13]. The smectic- A layer configurations 
corresponding with this splitting process is shown in Figure [6) 

The 25 layer simulation, with a more refined mesh, was used to both verify 
mesh independence of the 50 layer simulation (see Section 13.2.21) and obtain a 
more refined view of the smectic disclinations observed in this work. In order 
to identify disclination defects, the degree of biaxiality (see Section [^^ can be 
computed as follows [34,35]: 



f3' 



1-6 



[iQ-Q)-Q] 



(16) 



(Q : QT 

The core of the -1-1 disclination is shown in Figure[7^ as was determined from the 
25 layer simulation. The simulated structure of this disclination is both topo- 
logically correct and its biaxial halo/uniaxial center structure well-reproduced. 
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compared to its nematic counterpart [35-39]. The spherulite length scale of 
the 25 layer simulation was not large enough to observe the splitting of the +1 
disclination into a pair of +i thus Figure [7)3 shows the partially biaxial core 
structure of the +i computed in the 50 layer simulation. The core of this defect 
is also in agreement with its nematic counterpart, with a solely biaxial core [33]. 

3.2.2 Self-similar growth 

The radius of the spherulite versus time was determined for both simulations. 
Figure EK shows this data and a comparison of the results from the 50 layer 
and more refined 25 layer simulations show good agreement. Deviation in the 
early stages of growth result from the more refined 25 layer simulation capturing 
the formation and core structure of the +1 disclination more accurately. The 
convergence of the evolution of the two spherulite radii establishes the accuracy 
of the 50 layer simulation results beyond the radius of convergence. 

Following the shape dynamic growth regime, the spherulite forms a scale- 
independent, or self-similar, shape [32] . Spherulite growth obeys a power law r = 
i" in this growth regime, where n is both experimentally observed [40,41] and 
theoretically predicted [26,42] to approach n = 1 for phase-ordering processes 
quenched below the lower stability limit of the isotropic phase. A power law 
fit of the final two points of the spherulite radius evolution (Figure [H 50 layer 
simulation) shows that the simulation results in this work correctly predict a 
self-similar growth regime, where n — 0.95 is found. Additionally, the results 
from Figure[8]show that the transition from the initial shape dynamics regime to 
the self-similar growth regime occurs at a relatively small length scale compared 
to the growth of initially radial textured nematic spherulites [26] . 

Figure [9] shows the order parameter profiles during the beginning of the self- 
similar growth regime (Figure [5jii). Good agreement is observed between the 
bulk values of the order parameters (determined in Section [3. ip and those found 
in the dynamic simulation. The imposition of layer compression/expansion and 
curvature from the spherulite morphology and interfacial anchoring result in 
a decrease of the smectic-A order from that of the ideal homogeneous layer 
configuration. Additionally, a local increase in smectic-A order is observed 
in the region immediately outside of the spherulite core where layer compres- 
sion/expansion is minimized by the presence of the defect structure within the 
core (See Figure |6jii) . This phenomenon is an example of the competition be- 
tween gradients in long- and short-range order. The existence of such high 
smectic-A ordering in the presence of substantial layer curvature is in agree- 
ment with both theoretical and experimental observations that layer bending is 
a low energy distortion in the smectic-A mesophase [11,43]. 

4 Conclusions 

A comprehensive approach to modeling the first-order isotropic/smectic-A phase 
transition was presented and applied to phase diagram computation and growth 
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Figure 7: a) (left) t — 9.2/is (corresponding to the spherulites from Figures [SJ-ii, 
actual data taken from the 25 layer simulation) the biaxial halo and uniaxial 
center of a single +1 disclination observed b) (right) t = 2S.7^s (corresponding 
to the 50 layer simulation, Figure [Hii) the biaxial core of a +i disclination 
observed after splitting [14]; the surface corresponds to the degree of biaxiality 
(|16p , arrows correspond to the uniaxial director (and should be considered head- 
less), and the domain length scale is approximately / = Anm, or approximately 
one smectic layer, for both figures. 
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Figure 8: Log-log plot of the spherulite radius versus time for the 25 layer 
(dotted line) and 50 layer (solid line) simulations where a power law fit (r — t") 
of the final data points (50 layer simulation) yields n = 0.95. 
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Figure 9: t — 28.7^3 (corresponding to the 50 layer simulation, Figure [Hii) 
Order parameter profiles of S (solid line) and ip (dotted line); the symmetry 
axis perpendicular to that formed by the +-^ disclination pair was used and bulk 
values determined in Section [XT] for T = 330.0K are S = 0.75 and -0 = 0.20. 
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of an initially radial smectic-A spherulite in an isotropic phase. Summarizing 
the results determined from this work are as follows: 



• 



An optimized method of phase diagram computation was presented and 
applied (Section l2.4p . showing good agreement (Figure [3]) with the exper- 
imental data used in parameter determination (Section [ 



• Shape dynamics in the early spherulite growth period were found in agree- 
ment with the past theoretical work on a similar system, the first-order 
isotropic/nematic transition [14,24-28,30], where the initial imposed isotropic 
core forms a self-selected -1-1 disclination. As the spherulite radius in- 
creases, this -f 1 disclination splits into two -t-j disclinations, which is gov- 
erned by a competition of long- and short-range ordering in the presence 
of intcrfacial anchoring [13, 14] (Figures [5][7]). 

• A self-similar spherulite growth regime was observed following the ini- 
tial shape dynamic growth regime where a power law fit of the spherulite 
radius was shown to approach n = 1 (Figure[5]), in agreement with past ex- 
perimental [40,41] and theoretical [26,42] studies of mesophase spherulite 
growth under a "deep" quench. 

• The general structure of both -t-1 and -fi smectic disclinations were found 
in agreement with studies of their nematic counterparts [33,35-39]. Note 
that due to computational constraints, the mesh density used in this sim- 
ulation was adequate only to the extent of resolving the general structure 
of the defect core and not resolution with full detail. 

The simulation results presented show that the use of a high-order phenomeno- 
logical model and experimentally based model parameters results in a substan- 
tially more complete reproduction of the physical smectic-A system. Though 
current computational resources restrict this initial work to two-dimensions, 
these results show a strong correlation to past experimental [33, 35-41] and 
theoretical observations [14,24-28,30,42]. These promising results show that 
three-dimensional simulation of this model could be used to study the forma- 
tion and dynamics of smectic-A spherulites at length and time scales currently 
inaccessible via experimental study. 
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